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Abstract 

We provide a theory for the dynamics of cohapse of strongly charged poly- 
electrolytes (PEs) and flexible polyampholytes (PAs) using Langevin equa- 
tion. After the initial stage, in which counterions condense onto PE, the 
mechanism of approach to the globular state is similar for PE and PA. In 
both instances, metastable pearl-necklace structures form in characteristic 

4 

time scale that is proportional to A^s where A^ is the number of monomers. 
The late stage of collapse occurs by merger of clusters with the largest one 
growing at the expense of smaller ones (Lifshitz- Slyozov mechanism). The 
time scale for this process tcoll ~ A^- Simulations are used to support the 
proposed collapse mechanism for PA and PE. 
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Due to the interplay of several length scales charged polymers, referred to as polyelec- 
trolytes (PEs), exhibit diverse structural characteristics depending on the environment (pH, 
temperature, ionic strength etc) [|I|. The collapse of polyelectrolytes mediated by counte- 
rions is relevant in describing the folding of DNA and RNA. Reversible condensation of 
DNA induced by multivalent cations is required for its efficient packaging 0. Similarly, the 
role of divalent cations in enabling RNA to form its functionally competent state has been 
firmly established In both highly charged polyelectrolyte undergoes a collapse 

transition from an ensemble of extended conformations. 

Recently, simulations as well as theories have been advanced to describe the collapse 
of fiexible polyelectrolyte chains mediated by condensation of counterions An un- 

derstanding of the collapse transition in polyelectrolyte chains could be valuable in gaining 
insights into the condensation mechanism in biomolecules. Motivated by this we provide 
a description of the dynamics of collapse of fiexible PEs and polyampholytes (PAs). The 
theory for the latter becomes possible because of our suggestion (see below) that the physical 
processes governing the collapse of polyampholytes and polyelectrolytes are similar. 

In poor solvents (with reference to the uncharged polymer) the polyion is stretched at suf- 
ficiently high temperatures due to intramolecular electrostatic repulsion. The translational 
entropy of the counterions at high temperatures is dominant, and hence their interaction 
with the polyion is negligible. As the temperature is decreased the binding energy of the 
counterions to the polyelectrolyte exceeds the free energy gain due to translational entropy 
and they condense . This results in a great decrease in the effective charge on the polyion 
and, because the solvent is poor this leads to the compaction of the chain. 

These considerations have been used by Schiessel and Pincus to propose a diagram 
of states for highly charged polyelectrolytes. Consider a PE consisting of N monomers 
with a fraction 1) of monomers carrying a charge of — e. Strongly charged PE implies 
/(/b/^)^ > 1 where Ib = e'^ /AnekBT is the Bjerrum length, e is the dielectric constant of 
the solvent, and A is the mean distance between charges. At low temperatures we expect 
counterions (with valence z) to condense onto the PE because the Manning parameter 
(Z^M) > l/z. The electrostatic blob length C,ei is the scale at which the mutual repulsion 
between two charges is approximately /c^T, and is given by ^ei ~ Isz"^ /k"^ where k = —Incj) 
with (p being the volume fraction of free counterions [Q. The size of the PE is stretched i.e. 
L k'^A^N/lsz'^ provided ^ei < ~ ao0/(6 — T) where S,t is the thermal blob length, 9 
is the collapse temperature for the uncharged polymer. When < ^ei the PE undergoes a 
transition to a globular state. 

Here, we study the kinetics of this collapse process using Langevin equation. In our 
earher study we showed, using simulations, that the approach to the globular state occurs in 
three stages following a temperature quench 0]. In the initial stage the counterions condense. 
The intermediate time regime is characterized by the formation of metastable necklace-pearl 
structures In the final stage the pearls (domains) merge leading to the compact globular 
conformation. Here, the theory for the intermediate time regime is developed by adopting 
the procedure suggested by Pitard and Orland to describe the dynamics of collapse (and 
swelling) of homopolymers. We describe the dynamics of the late stages of collapse using 
an analogy to Lifshitz-Slyozov [|11[] growth mechanism. 

The equation of motion of the polyelectrolyte chain is assumed to be given by the 
Langevin equation 
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(1) 



where ( = kBT/D and D is diffusion constant of monomer, T is the temperature. The 
Hamiltonian of the ffexible polyelectrolyte chain with an effective Kuhn length is 

and 

Kw^r..r</.'^^<fM^ (3) 

Jo Jo |r(s, t) — r(s , t)| 

The thermal noise is assumed to be Gaussian with zero mean and the correlation is 

given by 

<7]{s,t)7]{s',t') >=2DS{t-t'). (4) 

In Eq.(0) r{s,t) is the location of monomer s at time t, Vc is Coulomb interaction between 
monomers at s and s', and Ib is the Bjerrum length. In writing Eq.(^) we have neglected 
hydrodynamic interactions p!0| , p!2| . Furthermore, excluded volume interactions are also omit- 



ted. The effects due to self-avoidance are not likely to be important in the early stages of 
collapse because in this time regime attractive interactions mediated by counterion conden- 
sation dominate. 

Collapse of Poly electrolytes: The polyelectrolyte chain is initially assumed to be in the 
G-solvent. At t = 0, we imagine a quench to a temperature below B so that the chain is 
effectively in a poor solvent. The equilibrium conformations of a weakly charged polyelec- 
trolyte (/(/s/A)^ -C 1), in poor solvents has been described by Dobrynin et al [|1^. and has 



been further illustrated by Micka et al g]. These studies showed that when the net charge on 
the polyelectrolyte chain exceeds a certain critical value then the equilibrium conformation 
resembles a pearl-necklace. This structure consists of clusters of charged droplets connected 
by strings under tension. This is valid when counterion condensation, which is responsible 
for inducing attraction between the monomers, is negligible. 

Consider a charged polyelectrolyte chain (/ = 1). Upon a quench to T < 9, such that 
the thermal blob length is not too small, the polyelectrolyte chain undergoes a sequence of 
structural changes enroute to the collapsed conformation [^]. In particular, after counterion 
condensation PE evolves towards metastable pearl-necklace structure. The dynamics of this 
process can be described using Eqs.(l-4) and the following physical picture. We assume that 
shortly following a quench to poor enough solvent conditions the counterions condense onto 
the polyelectrolyte chain. The time scale for condensation of the multivalent counterions 
is diffusion limited [Q, and is approximately given by tcond ~ p '^^^C/^bT where p is the 
density of the counterions. Explicit simulations reveal that the tcond is much shorter than 
time scales in which the macromolecule relaxes. [^]. Upon condensation of a multivalent 
cation with valence z{> 2), the effective charge around the monomer becomes {z — l)e. 
If the locations of the divalent cation occur randomly and if the correlations between the 
counterions are negligible then in the early stages, PE chain with condensed counterions may 



be mapped onto an evolving random polyampholyte |Tj]. We assume that the condensed 
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counterions are relatively immobile over the time interval tcond < t < tclust where tclust 
is the time required for the formation of pearl-necklace structures. With these assumptions 
the correlation between the renormalized charges on the polyanion can be written as 

mj = qlmj-'^)/iN-l) (5) 

where Qo is the net charge of the polyion after counterion condensation and the average is 
done over the conformations of the polyelectrolyte with condensed counterions. 

The above physical picture, which finds support in explicit numerical simulations [Q, can 
be used to calculate the dependence of the early stage dynamics of the collapse process on 
N using the method introduced by Pitard and Orland ||10|- The basis of this method is to 
construct a reference Gaussian Hamiltonian Hy{t) with an effective Kuhn length a{t) 

The value of a{t) is determined from the condition that the difference in the size of the 
polyelectrolyte chain computed using and the full theory vanish to first order in 5H = 
— H for all t. The equations of motion for the reference chain specified by r„(s,t) and 
for x(s, t) = r(s, t) — r„(s, t) are given (to first order in x) by 

where F{r{s,t)) = —d{Vc/kBT) /dr{s,t). The relaxation of each mode is obtained in terms 
of Fourier representation f„(t) = 1/N e'^'^'^^r{s,t)ds {un = ^) and is given by 

rn{t) = f dhGnit - ti)r/„(ti) + G„(t)f„(0) 

Xn{t)=D fdhGnit-h)[F4r,ih))-LuU^-^^)f^ih)] (8) 

where Gn{t) = exp{-D j' ^rft')- 

We assume that at t = the neutral chain is in G-solvent. Therefore < f„(0)fjli^(0) >= 
'^"'"i^mn- The thermal noise in Fourier space satisfies < fjnit) >= and < Vn{t)fj^(t') >= 



j^Snm^it — t'). We also assume that a{t) is a slowly varying function of time for t <^ trquse 
where trquse is defined implicitly by j^rouse _dr_ _ Thus, we can let Gn{t) = 
exp{—Duj'^t/ a?) for long wavelength modes. 

The time dependent R^g{t) can be expressed in terms of the amplitude and the relaxation 
rate of each mode 

Rlit) = E < kiwi >= E T^[(«o - a\t))G^{t) + a\t)]. (9) 

The variational parameter a{t) is obtained by equating the true value of the square of the 
radius of gyration to one calculated with H^. This implies that all the correction terms 
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to Rg computed using Eq.@ should vanish. To first order in xi^^t) the effective dynamic 



Kuhn length a{t) satisfies the self consistent equation [10| 







N 

<r^{s,t)x{s,t)>ds = (10) 



where the thermal average is done with respect to the noise term in the Langevin equation. 
The self consistent equation for a{t) is 



1_ _ ^ _ En r dt'Gnjt - t') < f^jt) > 



(11) 



o" 



where AH is the total interaction energy (Coulomb, excluded volume) of the PE. 

In the early stage of the coil-to-globule transition t <^ tclust, we can assume a^(t) ~ 
The initial driving force for collapse is the counterion-mediated coulomb attraction. All 
other forces (such as hydrophobic interactions due to the poor solvent quality) play a less 
significant role. Thus, we can set AH = Vc- With this observation the evaluation of 
< r„(t)F*(r„(t')) > gives 

< r„(t)F*(r„(0) > = < UtKit') > f f ^^W^^^^dsds' (12) 



37r2 (fi(s,s',t')/2)3/2 

with c„(s, s') = e*'^* — and Vl{s, s', t) = cf{s, s') < ff{t) >. The physical picture re- 
lating the state of the PE chain (till formation of pearl-necklace structure) to polyampholyte 
is used to perform an additional average over the "quenched" random charge variables (see 
Eq.(^). The resulting equation when substituted into the right hand side of Eq.(|llD yields 
an expression for the time dependent Kuhn length a{t), namely, 

a\t) = al[l - = a^l - (A)3/4) (13) 

Since qq' = —ql/{N — 1) when s 7^ s', we obtain Tc ~ {ao/lBY^^Cil/ D and it is nearly 
independent of the polymer size N . Note that is an estimate for the time scale in which 
the clusters in the metastable pearl-necklace structures form. Because the formation of such 
clusters occur locally (i.e. by interaction between monomers that are not separately by a long 



contours) Tc is expected to be independent of A^. Therefore, we identify Tc = tclust iH 
In the time range tcond < t < tclust, the radius of gyration is obtained by substituting 
Eq.(p!3D into Eq.(|]). We find the decay of Rg(t) can be approximated as 

Rl{t)^Rimi~{t/TPEr-) (14) 

with apE = 5/4 and tpe = i^^)^ (tclust)^ ■ The characteristic time tpe is therefore 

TpE ~ {Nao/lB)Hao/D)l. 

The formation of local clusters has also been suggested as a mechanism for collapse of ho- 
mopolymers in poor solvents. ]TB| , pT| The time scale for cluster formation for homopolymer 



has been calculated by Pitard and Orland [|TOl who found thp ~ {al/ D){al/\v2\Y^^N^/^ . The 
counterion-mediated attraction leading to pearl-necklace like structures occurs on a shorter 
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time scale tpe ^ N^/^ . In very poor solvents {Ib/O'o < |^2p/'2o)) where the hydrophobic in- 
teractions dominant the driving force for collapse of PE, we expect the dynamics to resemble 
that of homopolymer ||^ . 

Dynamics of PA collapse: The theory developed for PE is directly applicable to describe 
the dynamics of collapse of PA without having to justify the averaging over the quenched 
random charges (cf. Eq.(|^)). The Langevin dynamics for PA is given by Eqs.(l-4) where the 
quenched random charges explicitly satisfy Eq.(^). It is known that when the total charge 
of the PA Q < ±\/Ne the chain is collapsed at low temperature regardless of the quality of 
the solvent ||18|. The mathematical description of the dynamics given by Eq.(n) describes 



the early stages of collapse of PA. For PA, described by Eqs.(l-4), there is no counterion 
condensation. The metastable pearl-necklace structure forms after the random charges on 



the monomers are turned on provided Q < ±v A^e ||18[. We predict that the time for forming 
the metastable pearl- necklace conformations is also given by Eq. (|T^ The power law decrease 
in radius of gyration (Eq.(|l^)) is limited only to time scale within which the structure leads 
to several clusters. 

Late stages of collapse: For both PA and PE, the pearl-necklace structures merge (see 
below) at t > tqlust to form compact collapsed structures. This occurs by the largest 
cluster growing at the expense of smaller ones ("packman" effect) which is reminiscent 



of the Lifshitz-Slyozov |TT[ description of the kinetics of precipitation in supersaturated 
solutions. The driving force for this growth is the concentration gradient across the clusters 
or domains. If this analogy is correct then we expect that the size of the largest cluster S{t) 
to grow as 

S{t) ~ (15) 

with a |. The collapse is complete when S{t) ^ Rg{t oo) ~ aoN^ which implies that 
the characteristic collapse time tcoll ~ A^- 

In order to validate the proposed mechanism, i.e., the formation of necklace-globule 
structures and the growth of the largest domain by devouring the smaller ones, we per- 
formed Langevin simulations for strongly charged flexible PEs and PAs in sufficiently 
poor solvents so that the equilibrium structure is the compact globule. In Fig.(|I]) we display 
examples of the conformations that are sampled in the dynamics of approach to the globular 
state starting from 6-solvent conditions. Both panels (top is for PEs and the lower one is 
for PAs) show that in the later stages of collapse the largest clusters in the necklace-globule 
grow and the smaller ones evaporate. This lends support to the proposed Lifshitz-Slyozov 
mechanism. 

In order to estimate a (cf. Eq . ([T5|) ) we have used molecular dynamics simulations of 
polyelectrolyte and randomly charged polyampholyte and calculated the number of particles 
that belong to the largest cluster Nsit) as a function of time. In Fig.(||), the simulation 
results show the linear increase of Ns{t) for times greater than t ^ 30r for PA and t ~ lOOr 
for PE. Here r = al/D. Therefore, we estimate a ~ |. The change of slope for long times is 
due to the finite size effects and indicates the completion of the globule formation. Because 
the quality of solvent (see legend in Fig.(^ is not the same for PA and PE the onset time 
for linear behavior occurs at different times. 

Conclusions: We have presented a unified picture of collapse dynamics of PE and PA. 
Although the morphology of collapsed structures for PA and PE are different, our theory 
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shows that the mechanism of approach to the globular state for both should be similar. In 
particular, both PA and PE reach the collapsed conformations via metastable pearl-necklace 
structures. For PE the driving force for forming such structures is the counterion-mediated 
attractions. Charge fluctuations in PA lead to pearl-necklace structures. The life time of such 
structures is determined by a subtle competition between attractive interactions mediated by 
counterions and hydrophobic interactions determined by the solvent quality. The interplay 
between these forces may be assessed by estimating the free energy of necklace-globule 
structures. The necklace-globule conformation consists of n globules with nearly vanishing 
net charge that are in local equilibrium. The free energy of i*^ globule is 

All 

F,r.—{Af)R^ + 4naR^ (16) 

where Ri is the radius of the i^^ globule and a is the surface tension. Note that A/ is 
the same before and after the merger of clusters. The free energy difference between a 
conformation consisting of two clusters and the conformation in which they are merged is 
AF ~ (Sttct— 4(2)^/^7ro-)(A^/r;,)^/^. Typical charge fluctuation in each globule is ~ qo{N/ny^'^. 
If the coulomb energy fluctuation of each globule {SEfiuct ~ Qo/'^o{N/nY^^) is less than the 
free energy difference between the conformation with two separated clusters and the one 
where they are merged, then the system spontaneously grows to a large cluster. If the solvent 
quality is not very poor, i.e.~ (Svrcr — 4(2)^/^710") < ql/ao, the life time of the metastable 
necklace-globule can be long so that the collapse mechanism is controlled by the energy 
barrier between the metastable state and the globule. In this case, the attractive interaction 
between the clusters are induced by the mobile charge (counterion) fluctuations [jl9[ just as 



for the counterion mediated interaction between like charged rods. The morphology of the 
final collapsed state is also determined by a competition between electrostatic interactionis 
and hydrophobic forces. When the solvent is very poor the collapse state is amorphorous 
whereas when collapse is determined by AF the globular state is an ordered Wigner Crystal 
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FIGURES 




FIG. 1. Snapshots of the conformations sampled in the dynamics of approach to the globular 
state starting from S-solvcnt conditions. The final temperature is Ib = 5ao- For PA dark circles 
have q = +e and grey circles have q = —e. The simulations are performed by integrating the 
Langevin equations. The details are given Ref.[4]. The top panel is for PE and the lower one is 
for PA; N = 240 for both PE and PA; V2 = -2.97a^ for PE V2 = -6.18a^ for PA (the second 
virial coefficient V2 gives an estimate of the quality of the solvent [4]). The counterions are explicitly 
included only in the simulation of PE. The solvent quality is such that the combined systems of the 
collapsed PE chain and the counterions form a Wigner Crystal[4]. The corresponding morphology 
of PA is clearly more disordered (compare structure 'e' in the top and bottom panel). 
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FIG. 2. Time evolution of the number of particles belonging to the largest cluster Ns{t) in a 
randomly charged single polyampholyte and a polyelectrolyte (A^ = 240). The time scale for PA 
is given in the upper line. There are inherent numerical difficulties in computing Ns{t) for finite 
sized systems. For a given trajectory there is not a sharp boundary for the largest cluster because 
of this the error-bars are large. Nevertheless, we observe linear increase of Ns{t) for t > lOOr 
(PE) and t > 30t for (PA) where r = a^/D. This implies the Lifshitz-Slyozov growth mechanism 
i.e. S{t) ~ Ns{t)^^'^ ~ t^^^. At long times Ns{t) saturates indicating the completion of globule 
formation. 
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